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Abstract. The present work investigates the calculation of absorption and emission cyclotron line profiles in the 
non-relativistic and trans-relativistic regimes. We provide fits for the ten first harmonics with synthetic functions 
down to 10~* of the maximum fiux with an accuracy of 20 % at worst. The lines at a given particle energy are 
calculated from the integration of the Schott formula over the photon and the particle solid angles relative to 
the magnetic field direction. The method can easily be extended to a larger number of harmonics. We also derive 
spectral fits of thermal emission line plasmas at non-relativistic and trans-relativistic temperatures extending 
previous parameterisations. 
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1. Introduction 

The cyclotron-synchrotron radiation is produced by 
charged particles spiraling in a magnetic field. The non- 
relativistic limit (for particle speed /3 = w /c <C 1) is the cy- 
clotron radiation. The relativistic limit (for particle speed 
/3 = z;/c ~ 1) is the synchrotron radiation. At intermedi- 
ate speeds, the process is sometimes called gyro-magnetic 
(for /3 = z;/c ~ 0.5). 

The cyclotron-synchrotron effect is one of the most impor- 
tant processes in astrophysics and has been invoked in so- 
lar, neutron star physics and in galactic or extra-galactic 
jets. At non relativistic and ultra relativistic limits, the 
particle emissivity is provided by well known analytical 
formulae. 

The main contribution to the cyclotron radiation comes 
from low harmonics. In most astrophysical objects, this 
radiation is either self-absorbed or absorbed by the envi- 
ronment. It can even be suppressed by plasma effects such 
as in the Razin-Tsytovich effect. However, as stressed by 
different authors (Ghisellini & Svensson 1991, hereafter 
GS91; Mason 1992; GUozzi et al 1996) the cyclotron cross- 
section in the non- and mildly relativistic regimes is several 
orders of magnitude above the Thompson cross section. 
The cyclotron mechanism appears then as a very efficient 
thermalisation process in hot plasmas surrounding com- 
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pact objects. At higher energies, the synchrotron photons 
are a supplementary source for the Comptonisation pro- 
cess in accretion discs coronas, see e.g. Di Matteo et al 
(1997), GhiseUini, Haardt & Svensson (1998), hereafter 
GHS98, Wardzinski & Zdziarski (2000), hereafter WZOO 
and Wardzinski & Zdziarski (2001). 

A detailed numerical investigation of the thermalisation 
process in the presence of cyclo-synchrotron photons has 
scarcely been discussed (see however GHS98) but is of 
great importance for the spectral and temporal model- 
ing of the complex spectral energy distribution produced 
in compact objects. Low temperature plasmas that may 
be found in Gamma-ray bursts (in the comoving frame) 
and in accretion discs may produce cyclotron lines com- 
bined with synchrotron signatures. Such lines features can 
be produced in neutron star magnetospheres probing the 
magnetic field. 

The present work aims to propose useful fitting formulae 
for the first ten harmonics dominating the emissivity for 
particle velocities (or normalised momentum) between 
(3 = 0.01 and p = /37 = 2 (see figure 2 in Mahadevan, 
Narayan & Yi 1996 hereafter MNY96). We consider the 
case of an isotropic electron distribution in a tangled mag- 
netic field. The method can be extended to an arbitrary 
number of harmonics. At higher energies, however, the 
synchrotron formulae are worth being used. The derived 
functions can then be easily inserted into radiative trans- 
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fer codes. 

The article is organised as follows: Sect. [21 recalls the 
derivation of cyclo-synchrotron emissivity. Sect. |31 deals 
with the procedures used to compute the emission and 
the absorption coefficient: an integration over a small fre- 
quency interval around the cyclotron resonant frequency 
and a direct integration over the angle between the pho- 
ton and the magnetic field. These procedures are tested in 
Sect.0] Sect. El provides the fitting functions in the non- 
thermal case, at a given particle energy. Sect. provides 
the fitting functions for the emission coefficients produced 
in a thermal plasma. 



The ratio of two successive harmonics scales roughly as 

In the relativistic limit (/? — > 1) the emissivity is domi- 
nated by the harmonics of order n ~ 7^. The spectrum 
tends towards a continuum as the frequency interval be- 
tween two harmonics is Vb/l- The synchrotron power ra- 
diated by one particle per unit frequency interval is given 
by (we still assume 9p = 7r/2) 



dP 
dv 



K5/3(t) dt 



(5) 



1. 



x°-^ exp( 



-x) 



2. Cyclo-synchrotron emissivity 

The cyclo-synchrotron spectrum produced by one particle 
of mass m, charge e and velocity v = (3.c embedded in 
a uniform magnetic field B can be derived from classical 
electrodynamics (see for example Bekefi 1966, chapter 6). 
The emitted power spectrum, the so-called Schott formula 
is expressed as a sum of harmonics peaking at a frequency 
u = n VbM ('^ ^ l)j where 7 = (1 — is the Lorentz 
factor of the particle. The cyclotron frequency is Vb ~ 
e B/{2i: mc), where the magnetic field B is expressed in 
Gauss units. 



dP 
diydfl 



c 



where the critical frequency is Vc = 3/2 Vb^'^, x = vjVf. 
and K^ii,({) is the modified Bessel function of order 5/3. 
The particle emissivity scales as v^l'^ for x <C 0.29 (0.29 is 
the energy of maximum emissivity) and is exponentially 
cut-off at high energies (Ginzburg & Syrovatskii (1969)). 

3. Calculation procedures 

We developed two methods to numerically estimate the 
cyclotron emissivity and absorption cross-section of one 
particle at a given energy and frequency. Each method 
aims to treat the resonant term inside the Dirac function 
in a correct way, see Eqs. and (j^J. 
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b{y) , (1) Frequency integration 



with x — V 7/3_L sinO/vb, /3|| = f3cos6p, P± = PsinOp and 
Jn, J'n being the Bessel functions of (integer) order n and 
its first derivatives. 
The cyclotron resonance is defined by 



y — n— — 1/(1 — cos ( 
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The first procedure consists in integrating the expression 
(jH) over V in a, narrow frequency range [v — Aj//2, v + 
Az//2], with At/ — a u <^ V, a being a positive real vari- 
able. 

Eq. ^ becomes 
dP 



dv dfl 



(Op) 



a v c 
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Eq. Q is expressed in CGS units (e.g. ergs~^IIz~^st~^) 
as well as all equations in the text. 

The angles are defined relative to the dominant magnetic 
field direction: 6p is the particle pitch-angle and 9 is the 
photon angle. 

Integrating Eq.^over the photon angle and the frequency, 
averaging the pitch-angle leads to the total power emitted 
by one particle 



cos 9 — /3|| 
sin 9 



j2(x,) + /3ij;2(x,) 



• (6) 



Pf. 



4 B"^ 
-fjT c (7/3)2 

O O TT 



(3) 



where the Thomson cross section is ctt — 6.65 10^^^ cm^. 
At the non-relativistic limit, only the first harmonics con- 
tribute significantly to the emissivity. For the particular 
pitch-angle 9p = 7r/2, and after integration over all ob- 
server angles, the power emitted in a given harmonic n is 
(Bekefi (1966)) 



The resonant term is Xr = n(3± sin6'/(l — cos6'). The 
limits ni and 112 are fixed by the resonant condition Q 
taken dX v — Ai//2 and v -\- Ai//2. 

3.2. Direct integration 

The second procedure eliminates the resonance through 
a direct integration over 9, the angle between the photon 
and the magnetic field. The resonance condition selects an 
angle for a given harmonic n, a given pitch-angle and (3 



cos 9r 



71/ 



(7) 



_ Stt^ e'vl (n + l)(n2"+i) 
^" " c (2n+l)! ^ 



(4) 



Integrated over the solid angle, Eq. leads to 
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The resonant term is now Xr = sind^jv^,. The sum 

boundaries n\ and ni define the range of n for which 
|cos6lr|<l in Eq.d 

We note that for cos 6*^=0, cos6'r becomes undefinite 
and Eq. © breaks down. In this case, all the power is radi- 
ated at the resonant frequencies, i.e. dP/dv{9p = 7r/2) = 
except at frequencies v — mvf)/^ and at those frequencies 
the power per unit frequency is actually infinite. The to- 
tal power emitted in harmonic n at pitch angle Op — -k /2 
may be estimated by integrating Eq. over frequencies 
around the resonance (leading to a result similar to Eq. 
and then numerically over the observer angles 9. 



3.3. Discussion: total line emission 

The power spectrum dP/dvln at a given harmonic is ob- 
tained by integrating numerically Eq. (O over the photon 
solid angle and over the particle pitch-angle ftp using 
an extended Simpson formula. For the direct integration 
in Eq. (jSj) only one integral over the particle pitch-angle 
is necessary. 

Both methods fail to account for the unavoidable reso- 
nance at frequencies ~ ni^t/j that is due to the con- 
tribution from pitch angle 6p = tt/2. In the frequency in- 
tegration method, the Dirac function is artificially broad- 
ened over the integration bin of width av. This leads to 
an overestimate the line emissivity around the resonances, 
but the resulting frequency integrated power is conserved. 
In the direct integration scheme we neglected the reso- 
nant contribution at frequencies v = nvb/"/- Although this 
procedure gives more accurate line profiles, the frequency 
integrated power is slightly underestimated. The two in- 
dependent methods are compared for different particle en- 
ergy regimes in Fig. ^ They appear to be in agreement 
within ~ 10% at low energies and ^ 20% at high par- 
ticle energies. The errors (the relative difference between 
the two emissivities) have been calculated down to 10^^ 
the line maximum. The calculated flux obtained from the 
second method appears to be lower around the peak of 
emission at least for the first few harmonics and at low 
particle energies. The discrepancy is always of the order 
of 10%, less than the accuracy we aim to reach using the 
fitting functions. For this reason, we will only consider the 
first (frequency integration) method for the tests detailed 
in Sect. 13 For both methods, the numerical integration of 
the power spectrum over frequencies gives the total power 
(Eq. O within a few percent. 

The best agreements (between the two methods) were 
obtained using an a coefficient (associated with the fre- 
quency integration method) in the range 10^'^ to a few 
times 10~^. These values ensure a stable numerical inte- 
gration. MNY96 have developed a similar but different 
method, approximating the Dirac peak using a smooth 
function over a frequency interval av. They reach a good 
numerical stability with a = 5 10~^. In the following, we 
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Fig. 1. Relative differences between the two formulae (O 
and deriving the power emitted by one particle. The 
energy is scaled to the line peak energy. The harmonics 1 
and 5 are shown at three particle energies (3 = 0.1, 0.5, 0.9 
corresponding to the energy ranges explored in this work 



use the first method to find the fitting functions at a given 
particle energy (see Sect. E)), the second method has the 
advantage of avoiding one integration and will be used in 
the thermal plasma case (see Sect-EJ. 

3.4. The cyclo-synchrotron absorption cross section 

The cyclo-synchrotron absorption cross section has been 
derived by GS91 from the Einstein coefficients in a three 
level system. The differential cross section is the differ- 
ence between true absorption and stimulated emission 
cross sections. In the relativistic limit, the stimulated ef- 
fect takes over for large photon angles with the magnetic 
field, but the total cross section is always positive. 
If the photon energy is much smaller than the electron 
kinetic energy ^ ^m^c^ / h) , the angle integrated cross 
section can be expressed by the angle averaged emissivity 
Jii^jP) — {dP/dv)/ATT in a differential form as 



1 



1 d[-fP3{u,p)] 



2 rriei' 7 p 



(9) 
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The cyclo-synchrotron cross section can reach values ~ 
lO^^ar/BG for u = i^b, and P ~* 0. 

The absorption cross section at a given frequency for a 
given harmonic (riQ < 10, f3 and p are fixed) is calculated 
using Eq. 0. While integrating over the solid angles, we 
first test if the harmonic uq is in the range of the permitted 
harmonics evaluated at a given pair of angles 0,9p. The 
harmonic number interval is defined from Eq.lPJ by 



m = E[{1 - a) 1/ 7 (1 - /3||Cos(6i))] 
n2 = E[{l + a) 1/7 (1 - /3||cos(0))] 



1 



(10) 



Once the test is successfully passed, we calculate the ana- 
lytical derivative over 7 of the emissivity at the harmonic 
number no, e.g. the derivative of 



dP 



2TTe'^vl 



a V c 7^(1 — /3|| cos6')^ 



is6'-/3|| 
sin 6 



JnoM+PljZM 



,(11) 



with Xr = nof3± sin 9/ (1— /3|| cos 6). We finally sum at each 
frequency bin the contribution from the ten first harmon- 
ics to get a{i^, P). 
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4. Tests 

We performed a series of tests for our numerical schemes 
in both mono-energetic and thermal cases. The methods 
developed are not well adapted to the relativistic regime 
as the number of harmonics to be summed increases rad- 
ically when [3^1. They failed to converge to the correct 
analytical expression Q at high frequencies. 
We first tried to recover the known cyclotron limit, the 
total power radiated by one particle, at a given harmonic 
and compared it with expression Q . The second test com- 
pares the numerical cyclotron absorption coefficient at a 
given harmonic with the non-relativistic coefficient pro- 
vided by GS91. For the thermal plasma case, we recover 
the cyclotron absorption coefficients given by Chanmugam 
et al (1989) (hereafter C89) and the Kirchoff law. 



Fig. 2. Relative differences between Eq. Q the power 
emitted by one particle integrated over the solid angle 
and Eq. Q the power radiated by one particle in the non- 
relativistic regime. Both emissivities have been calculated 
at a pitch-angle Op — tt/2. The plot is provided for the 
ten first harmonics and for two different particle energies 
(3 = 10^^ (triangles), and 0.1 (crosses) . The emissivities 
have been normalised to the harmonic one. For scaling rea- 
sons, the errors have been divided by 50 for the (3 = 0.1 
case. 



agreement and a cross section of the order of C7{v, (3 
0) - lO^^Bg^ (7t- In the trans-relativistic regime, we also 
found a good agreement at frequencies v > (see Fig. |S1 
lower panel). 



4.1. Tests for mono-energetic distribution 
4.1.1. Cyclotron limit 

In the cyclotron limit {/3 ^ 0), we integrate the first har- 
monics over the resonant frequency range with a particle 
pitch-angle 9p — tt/2 and compare our results with Eq. I^J. 
We plot the relative differences for different /3 for the first 
ten harmonics in Fig. 13 As can be seen, the agreement is 
quite good, even if the errors increase with the harmonic 
number, however they are limited to few percent. 
Concerning the absorption cross section we compared our 
results with the analytical expressions given by GS91. We 
point out that their formulae have to be corrected by a 
factor 1/4 TT in order to verify Eq. Then, at low en- 
ergy and at frequencies ~ Vb, we did obtained a good 



4.2. Thermal absorption coefficients 

C89 (and references therein) derived the thermal ab- 
sorption coefficients for non relativistic thermal plasmas 
(T < 50 keV) for the two polarisation modes. The 
ordinary mode has an electric vector field oriented 
along the magnetic field given by the first term in 
Eq. 1^, while the extraordinary mode has an electric 
vector field perpendicular to the magnetic field given 
by the second term in Eq. Q. The extraordinary mode 
usually dominates. At different harmonics, we calculate 
the absorption coefficients from the emissivity at a 
given angle 9 integrated over a Maxwellian distribution 
iVe/(e/-i:2(l/e)) 7^/3 exp(-7/e). We have defined the 
dimensionless temperature as = fc^T/m and A^e the 
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Fig. 3. Thermal absorption coefficients (in cm~^) at 
T = 5, 15, 25 keV. The photon angles are 6 — 
80° three dotted — dashed hne , 60° long — dashed hne, 
30° dot — dashed line. The results can be compared with 
C89. The calculations have been performed for a magnetic 
field B = 10 MG. 



electron density. 

We plot our results in figure 13 for three temperatures at 
three different 9 values and at different frequencies cov- 
ering the first 20 harmonics. The agreement with C89 is 
good (see their tables 6A, 6C and 6E), the numerical ther- 
mal coefficient falls close to the extraordinary mode val- 
ues. 

Using the same calculation procedure, we check that both 
thermal emission and absorption coefficients fulfill the 
Kirchoff law, ji, = Bi,{T) a^- We plot, in figure^ the 
relative differences between / ctu obtained from numeri- 
cal integration with B,^{T) versus the frequency for non- 
relativistic temperatures. The agreement is again good ex- 
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Fig. 4. Relative differences between the ratio of the 
calculated emission and absorption coefficients jv/oii, 
and B^{T) at non-relativistic temperatures T = 
1 (triangles), 10 (squares), 50 (stars) keV, in the frequency 
range 1-100 I'/i't- 



cept for T > 10 keV at ~ where the flux is weak and 
the uncertainties obtained with the numerical integration 
method increase. 



5. Non-thermal lines 

We now describe the fitting function for the mono- 
energetic case in the non- and trans-relativistic regimes. 
Whenever it was possible we used simple (essentially poly- 
nomial) functions defined by the smallest number of pa- 
rameters. Accounting for the evolution of the line profile 
with the particle energy, we split our parameters into two 
classes: the parameters below the energy of the line peak, 
the parameters above the energy of the line peak. The line 
normalisation with /3 and n is also parametrised. 
The generic form for the fitting functions is 



F{n,f3,iy) = Ni 



Ki-\ 



where 



'Anf 



- V 



for V < I/max , 
for V > i^max . 



(12) 

(13) 
(14) 



^ sup 



The different frequencies are = Vh njipi [1 -I- /?]), 

i^sup = i^b n/^'j [1 — /?]) given by the resonant condition 
Eq. The maximum frequency is Vmax — n/^. This 
choice respects the normalisation F{n, (3,Vmax) = Ni. 
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The parameters are chosen in order for the functions 
(|12() to fit the fine (where the worst error was 20 %) in a 
frequency range Vmax up (down) to a frequency i/gnd > 
{i^end <) where the hne flux is 10~^ the peak flux, we have 
the ordering i/i^f < {v^nd <) < Vmax < {>^end >) < v sup- 
Tins value was chosen to ensure good flux estimates 
at frequencies between two lines maxima where in the 
trans-relativistic regime, two subsequent harmonics may 
both contribute significantly to the emissivity. 
The parameters are A^i (n, (3) > Q for the line normalisa- 
tion, Ki{n,l3) > 0, /i(n,/3) > and hin.P) > for the 
line profile. 

Once a line is fitted, the n and (3 dependencies of the four 
primary parameters Ni,Ki,Ii, I2 have to be determined 
using other fitting functions. The whole fitting procedure 
is twofold, e.g. we first derive the previous coefficients in 
a grid {P,n) using Eq. ((T^ . We called this the primary 
fitting procedure. Wc then provide new fitting func- 
tions for both (3 and n dependencies of the four previous 
parameters, this is the secondary fitting procedure. 
Obviously, the general fitting function (|12|l is not 
unique, other forms may be much more accurate for 
the primary fitting procedure but they usually require 
more parameters and lead to a complex second step. 
The solutions obtained are polynomials of degrees 
lower than 4, they are accurate enough to be used in 
a wide range of problems dealing with cyclotron radiation. 

We consider two energy regimes: the cyclotron-gyro- 
magnetic (non-relativistic) regime defined by /? < 0.5 
and in the trans-relativistic regime defined by p < 2 and 
P > 0.5. Beyond this, the usual rclativistic formulae apply. 
This choice was motivated by the different behavior of the 
primary parameters in the two regimes (see next). 
For both emission and absorption lines, the second expo- 
nent, at a given harmonic, l2{n) has been fixed at /3 = 0.01 
(the lowest energy considered here) and has an imposed 
evolution in the non rclativistic regime (up to P = 0.5). 
Above, in the trans rclativistic regime, we re-initialise 
/2(n) to its value ai (3 = 0.5 and impose an other evo- 
lution with p. This procedure limits the fluctuations in 
the other two primary parameters Ii and Ki. We also 
had to consider the evolution of the first harmonic indi- 
vidually for both emission and absorption coefficients, and 
the evolution of the second harmonic also individually for 
the absorption coefficient. 

The number of parameters being high (6 parameters for 
the shape and one for the normalisation), the results are 
unfortunately too complex to be presented in tabular 
forms. We have therefore decided to provide fortran rou- 
tines calculating the power spectrum for frequencies i' < 
10 fb- This procedure has been adopted for the thermal 
case too (see next section). All the programs are accessible 
via ftp anonymous at ftp.cesr.fr/pub/synchrotron (in 
cyclosynchro.tar.gz). We simply present here results in a 
graphical form. 



5.1. Results 

For each figure, the first plot presents the numerical cal- 
culation and their fits superimposed and the second plot 
presents the relative errors. Figures El and El display the 
results for both emission and absorption lines in the non- 
relativistic limits at particle energies corresponding to 
P — 0.1 and P = 0.3. Figures and jSl display the results 
from the trans-relativistic regime at particle energies cor- 
responding to P = 0.7 and P = 0.9. 

The units are cm~^ and erg s~^ Hz"^ cm^'^ st^^ for 
the absorption and emission coefficients respectively. 

In most of the cases, the harmonics are fitted down 
to 10"'* the peak line flux with errors never exceeding 20 
% (ratio larger than 1.25 or lower than 0.8). Exceptions 
may be found in figures and for frequencies v > v^ax 
at harmonic number n > 4 where the profile is very 
sharp. The fiux is overestimated here but for a restrained 
frequency range. However, even in these cases errors are 
around 20 % and the higher harmonics do not strongly 
contribute to the overall fiux. The errors increase at the 
edges of a given harmonic in the non-relativistic regime 
and at frequencies v of the order Vinf or Vsup in the 
trans-relativistic regime due to the low accuracy. The 
noise appearing at the edge of the frequency domain 
(1-10 Vh) is due to the numerical integration method. As 
the fiux decreases, the method would have required a 
frequency dependent (increasing) a in Eq. ^ to keep a 
smooth line profile. As the fiux at these energies are far 
beyond lO"'* the peak line flux, we disregarded the prob- 
lem. However the emission (or absorption) beyond 10 Vb 
is not well reproduced since at these energies harmonics 
of higher order contribute and the synchrotron formula 
is worth to be used. Finally note that the flux in the 
absorption of the first harmonic in the trans-relativistic 
case is under-estimated by a factor two under Vmax for 
fiuxes lower than 10~^ the peak fiux. 
Also presented are the non-relativistic and the rcla- 
tivistic limits derived by GHS98 and GS91 respectively. 
In figure (HI one can see a good agreement at frequen- 
cies V > Vb between the numerical and the analytical 
calculations, for both emission and absorption coefficients. 



6. Thermal emission lines 

The previous emission coefficients are now integrated over 
a rclativistic Maxwellian distribution at a given temper- 
ature O = fcs T/mgC^. As stated above (Sect. j^J, for 
computational facilities, we used the second method (di- 
rect integration). 

We computed the line spectra for the ten first harmonics 
at 200 different temperatures forming a logarithmic grid 
in the range of O = 10^^-1. We make available electron- 
ically a table of emissivities in the 10 first harmonics as a 
function of photon frequency and a fortran routine read- 
ing this table and interpolating between the precomputed 
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Fig. 5. Results for non-thermal lines at jS = 0.1. The 
first upper panel shows the spectra produced using the 
fitting functions (solid line). The flux obtained from nu- 
merical integration indistinguishable at this scale has been 
skipped for a better leasibility. The ratio of the fitting 
function to the numerical results is plotted in the lower 
panel. The second panel shows the fits and the ratio 
for the absorption line. The emission and the absorp- 
tion coefficients arc normalised to {2TrneeB^) / {irieC^) and 
3/8 7r A^e {<^T/ctf) (i?cr/S) respectively (a/ — f/f37isthe 
fine structure constant and Bcr ~ 4.4 10^'^ Gauss is the 
critical magnetic field). The same normalisation is adopted 
for all figures. The formula (8) of GHS98 and (2.23) of 
GS9f in the non-relativistic limit has been used to plot 
the emission and absorption coefficients in dotted line. 



Fig. 6. Same as in figure El but for /? = 0.3. The formula 
(2.16) of GS91 in the relativistic limit for the emission 
coefficient is plotted in dotted line in the upper panel. 
The relativistic formula (2.17) of GS91 in the relativistic 
limit for the absorption cross-section has been used to 
derive the absorption coefficient plotted in dotted line in 
the lower panel. 

and 1^2 so that: 



emissivities. The absorption coefficients can then be de- 
rived from the Kirchoff law. 

We found that the emission due to an individual harmonic 
can be conveniently approximated as follows: 
Let Vmax be the photon frequency of the peak of harmonic 
n (see previous paragraph). We define the frequencies vi 



' {n + e^b 

V2 = (n + 0.001)1/6, 

we have < Vmax < V2- 

We introduce the following notations: 

X2 = 1^2/ {nvb)- 



(15) 
(16) 



(17) 
(18) 
(19) 



Let A^i, and Nmax be the line emissivity at vi, V2 and 
Vmax- We further note pi, p2, P3... etc, the parameters 
depending on n and Q that provide the best fit to our 
numerical results. Then, for y<vi., the thermal emission 
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Fig. 7. Same as in figureElbut for f3 — 0.7. 

coefficient can be conveniently approximated by the func- 
tion: 



Fig. 8. Same as in figure |S1 but for (3 — 0.9. 



f{ti, e, I/) = iVi 

exp 



1 -P2X1 
Vx J 1 — 'P2X 



For V2<v: 



\V2 



-P6 



+ X — xl 



(20) 



exp 



P8 

29 



X2 ~ I + prvQ 

- X~ X2 



~ 1 
1 

X2 



(27) 



29 \x Xl 

For yi<U<l',nax- 

F{n, 9, ly) = iVmaa; exp [a ln(t//t/,„aa;)^ + blniv/vrnax)^] (21) 
where : 

a = (3-p4)ln(iVi/iV„,a:r)/ln(!^i/i^™ax)^ (22) 

6 = (P4 - 2) ln(iVi/Ar„,a:r)/ln(z^l/l^™aa;)^ (23) 
For J^mai,<I^<I^2: 

F{n,Q,iy) = N^axexp [c\n{v/vmaxf + dlii{iy/iymaxf] (24) 
where : 



C = (3-p5)ln(iV2/A'"maa:)/ln(z^2/!^maa;)^ 
d = (P5 - 2) ln(A^2/A^ma:r)/ln(z/2/j^ma:r)^ 



(25) 

(26) 



This fitting function thus depends on 12 parameters 
(i^max, Nmax, Ni, N2 plus 8 parameters Px) that were de- 
termined from fits to our tabulated emissivities. We then 
constructed a table giving the best fit parameters as a 
function of temperature and harmonic numbers. Contrary 
to the case of the non-thermal emission, we have not been 
able to find simple fitting functions approximating con- 
veniently the temperature dependence of these parame- 
ters (the second fitting procedure described in Sect. 
However, the table of parameters we derived may be used 
to directly estimate the thermal synchrotron emission and 
absorption coefficients at low photon frequencies. We pro- 
vide electronically the parameter table together with a 
fortran routine computing the coefficient from that ta- 
ble. This methods is computationally faster and much less 
memory consuming than directly interpolating between 
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Fig. 9. From top to bottom: Thermal synchrotron emis- pig. 10. Thermal synchrotron emission coefficient for 
sion coefficients for kTe=5M IQ-^, 0.153 and 0.511 keV. 5.11 and 15.3 keV from top to bottom respec- 

Upper panels: the solid curves show the approximation us- tively. See caption of Fig. |21 
ing the fitting functions described in Sect. El The results 
from numerical integration are undiscernable. The dotted 
curves show the results from the formula of WZOO. Lower 
panels: ratio of the fitting functions to the numerical re- 
sults. 
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Fig. 11. Thermal synchrotron emission coefficient for 
fcTe^Sl.l, 153 and 511 keV from top to bottom respec- 
tively. See caption of Fig. El 

the tabulated emissivities. The solutions presented 
here are much more accurate than the parametrised func- 
tions obtained by different authors (MNY96, WZOO and 
references therein) at low temperatures and photon fre- 



quencies for which the errors around v\, are particularly 
severe. 

6.1. Results 

Figs. El El and m compare the emission coefficient in 
erg Hz~^ cm~^ st~^ obtained using the fitting func- 
tions with the results of the numerical integration at differ- 
ent temperatures. The approximation appears to be good 
for the whole range of temperatures. For illustration we 
also show in these figures the results of the analytic ap- 
proximation given by WZOO. As long as Q< 1 this analyt- 
ical estimate connects smoothly with our results around 
u ^ lOvb and may be used to compute the emission at 
higher frequencies. On the other hand, for 0^1 and 
higher, the approximation of WZOO is not so good, and our 
approach based on individual harmonics becomes inaccu- 
rate around v ~ lOvb due to the significant contribution 
from higher orders. 

7. Conclusion 

In this work, we presented two independent methods to 
compute the cyclo-synchrotron emission and absorption 
coefficients in a tangled magnetic field in the non- and 
trans relativistic regime. These calculations were done 
both for mono-energetic and thermal electrons. Only the 
first ten harmonics dominating the flux in theses regimes 
have been considered, but the procedure can easily be ex- 
tended to a larger number of harmonics. We found syn- 
thetic fitting functions that reproduce our numerical re- 
sults with a good accuracy (with a relative error at worst 
20%). They complement the previously existing approx- 
imations valid only at higher photon and electron ener- 
gies. Although these fitting functions are probably too 
complicated to be used in analytic calculations, they may 
prove very useful for numerical purposes. For instance, 
they could be used in high energy radiation transfer codes 
including both photon and particle dynamics (e.g. Coppi 
1992; Stern et al. 1995). The fitting functions are made 
available electronically in the form of tables and fortran 
routines that can be easily implemented in such codes. 
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